Rotational quenching of CO due to H2 collisions 
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ABSTRACT 

Rate coefficients for state-to-state rotational transitions in CO induced by 
both para- and ortho-H2 collisions are presented. The results were obtained us- 
ing the close-coupling method and the coupled-states approximation, with the 
CO-H2 interaction potential of Jankowski & Szalewicz (2005). Rate coefficients 
are presented for temperatures between 1 and 3000 K, and for CO{v = 0,j) 
quenching from j' = 1 — 40 to all lower / levels. Comparisons with previous 
calculations using an earlier potential show some discrepancies, especially at low 
temperatures and for rotational transitions involving large \Aj\. The differences 
in the well depths of the van der Waals interactions and the anisotropy of the two 
potential surfaces lead to different resonance structures in the energy dependence 
of the cross sections which influence the low temperature rate coefficients. Ap- 
plications to far infrared observations of astrophysical environments are briefly 
discussed. 

Subject headings: molecular processes — molecular data — ISM: molecules 
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INTRODUCTION 



Molecular hydrogen and carbon monoxide are the most abundant molecular species in 
the majority of interstellar environment s, and they p lay important roles in determining the 



physics and chemistry of diffuse clouds (IHerbst 



200 ll ). CO is observed in diffuse clouds in 



absorption in the ultra violet (UV) and visible and in emission at infrared (IR) wavelengths 



(ISnow fc McCall 



20061) ■ revealing t 
sight. As an example, iLisztl ( l2006l . 



le ph ysical and chemical complexity along these lines of 
20071 ) recently investigated the formation, fractionation, 
and rotational excitation of CO in H2 bearing diffuse H I clouds. Due to its small rotational 
constant, CO can be easily rotationally excited by collisions with other species in interstellar 
gas, mostly H2, with CO rotational lines providing important diagnostics of gas density and 
temperature. The abundance ratio, CO/H2, is assumed to be roughly constant in dense 
molecular gas, as such, CO is often used as a tracer of H2, as the latter is difficult to detect 
in emission in cold environments. The observed CO abundance is therefore used to estimate 



the total H2 content fiSonnentrucker et al. 



20071). 



In environments with an intense UV field, the radiation can drive the chemistry and 
internal level populations out of equilibrium. In such situations, a photodissociation region 
(PDR) resides at the interface of the hot H II region and the cold molecular region. In recent 



with an emphasis 


on the excitation mechanisms of both 


iy, anc 


I CO ( 


van Dishoeck & Black 


1988; 


Warin et al. 


1996; 


Le Petit et al. 


2006; 


Shaw et al. 


2005 




RoUie et al. 


2007 


)■ 



In more dense environments such as low-mass dwarf stars (i.e., M, L, or T dwarfs), 
CO is an i mportant opacity source and its infrared line list can be obtained quite 



accurately (jPavlenko fc Jones 



2002 



Jones et al. 



2OO5I ). It has also been observed in the 



dayside spectrum 



( ISwain et al. 



of the transiting hot- Jupiter extrasolar giant planet (EGP) HD 189733b 



20091). While molecular level populations are typically treated in local 
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thermo dynamic equilibrium (LTE) in stellar and planetary atmospheres, ISchweitzer et al. 



( 120001 ) have pointed out that non-LTE conditions may exist in very cool stellar atmospheres 
or in cases where the temperature of the incident radiation, i.e. close-in EGPs, is much 
different from the kinetic temperature of the gas. 

Due to the astrophysical importance of molecular hydrogen and ca rbon monoxide, the 



CO-H^ col 



Butz et al 



isiona l system has been t 



he su 



1971 



Kuppermaii et al 



1973 



ect of numerous experimenta . 



Bouanich fc Brodbeck 



1973 



Kudian et al 



1967 



Nerf &: Sonnenberg 



19751: iBrechignac et al.l Il980l: lAndresen et al.l Il982l: 



icard-Bersellini et al. 1983: 



2000 ) and theoretical ( 



Poulsen 



1982 



Schramm et al. 1991; Drascher et al. 19981: McKellai 1990. 1991, 1998; Antonova et al. 



Green fc 



Schinke et al 



1984; 



laddeud 1976 



Bade et al. 



Prisette et al. 



1985a b 



1978 



Flower et al. 



Jankowski fc Szalewicz 



1998, 



979 



2005; 



van Hemertlll983l: iParish et al.lll992l: iDanbv et al.lll993l: ISalazar et al.lll995l: iReid et al. 



1997 



Mengel et al 



Antonova et al. 



2001 



2000 



Baker fc 



Wernli et al 



'lower 



2006 



1984 



Flower fc Launav 



Yang et al. 



2006 



1985 



Flower 



2001 



b|) studies. Quantitative 



determinations of state-to-state cross sections and rate coefficients for CO-H2 collisions 
are crucial to numerical models of various astrophysical environments, such as those 
highlighted above. However, as measurements of these quantities are difficult, numerical 
models often rely on cr oss sections and rat e coeffi cients derived from theoretical methods. 



In a highly-cited work. 



Green fc Thaddeud (Il976[ ) performed close-coupling calculations of 



rate coefficients based on an approximate CO-H2 potential surface. Since then, a number 
of quantum scattering calculations were carried out on various potential energy surfaces 
(PESs). In this work, we have performed comprehensive calculations of state-to-state cross 
sections in an effort to develop a complete database of rotational quenching rate coefficients, 
but utilizing a more recent, and presumably more accurate PES. In the following sections, 
the choice of the potential surface, the adopted scattering approach, and the results are 
discussed. We briefly highlight important astrophysical applications of the current CO-II2 
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rate coefficients for modeling of far IR (FIR) observations. 



THE POTENTIAL ENERGY SURFACE 



1984 



1 



Bacic et al.Ml985a]l : jjankowski fc SzalewiczlilQQSl 



A number of PESs f^Prisette et al. 1978 : Flower et al. 1979t Poulseru 1982 : Schinke et al. 
^1 u ajj — ^ . u k 

20051 ) have been developed for the 



CO-H ? complex. One of the most accurate surfaces was given by 



Jankowski fc Szalewicz 



( 119981 ) who calculated a four-dimensional, rigid- rotor PES referred to as Vgg. This 



surface was constructed using the symmetry-adapted perturbation theory (SAPT) method 
with high-level electron correlation effects. An analytical fit of the ab initio PES has a 



global minimum of 109.3 cm ^ at the intermolec ular separation o: 



Antonova et al. 



7.76 ao for the linear 



teoOOh adopted Vgs in 



geometry with the C atom pointing toward H2. 
their calculations of state-to-state cross sections for CO rotational excitation at collision 
energies between 795 and 991 cm~^ which were found to be in good agreement with their 



measurements. The Vgg PES was also utilized in seconc 



mixtures of hydrogen and CO (IGottfried Sz McBane 



virial coefficient calculations for 



2OOOI). A comparison with experimental 



data suggested that the van der Waals well of the Vgg PES is too deep by 4%-9%, though 
it represents an improvement over earlier surfaces. The Vgg surface has also been used in 



full coupled-channel cross section and rat e coefficient c a. 



scattering of CO by para- and ortho-H2 (iFlower 



2001 



culations for rotationally inelastic 



Mengel et al 



2001). However , 



Mengel et al. 



(120011 ) 



because the attractive well of Vgg was presumed to be too deep, 
modified the PES by merely multiplying the surface by a constant factor of 0.93 which was 
subsequently used in their scattering calculations. 



A newer CO-H2 PES was reported by 



Jankowski fc Szalewica (120051 ) . referred to as 



V04. To improve the surface accuracy, they adopted the coupled-cluster method with single, 
double, and noniterative triple excitations [CCSD(T)] and the supermolecular approach. 
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The Vo4 surface was obtained on a five-dimensional grid including the dependence on the 
H-H separation, but with the CO molecule taken to be rigid with the C-0 separation set to 
the value of the C-0 distance averaged over the CO ground state vibrational wave function. 
The PES was then obtained by averaging over the intramolecular vibration of H2 to yield 
a four- dimensional rigid-rotor surface. Similar to Vgs, an analytical fit of the V04 was 



presented which has a 



obal minimum of 93.05 cm ^, or 16.22 cm ^ shallower than Vgg. 



McKellai] (119981 ) for para-Hs 



Jankowski &: Szalewica (120051 ) used V04 to calculate the rovibrational energy levels for the 
para- and ortho-H2-CO complex as well as the second virial coefficient . The rovibrational 
energies were found to agree with the experimental values of 
to within 0.1 cm^^ while a scaling of V04, resulting in a deepening of the well depth by 4 
cm~^, was needed to match the measured second virial coefficients. V04 may be too shallow 
though it is possible that the experirn ental second virial coefficients are systematically 



too low (Jankowski fc Szalewicz 



2OO5I ). Comparative studies of coUisional cross se ctions 



using the two surfa ces h ave been pr e sented in our earlier work (lYang et al 



Wernli et al. 



(12006h . In 



Yang et al. 



2006af ) and by 



(j2006al ). calculations with the V04 s urface gave cross 



sections in excellent agreement with the state-to-state measurements of 



Antonova et al. 



(l2000h . 



3. THEORETICAL APPROACH 



he t heory for the scat tering of two linear rigid rotors has been developed in 



(Il975h and 



HeiletaL 



Green 



(119781 ). The calculations presented here were performed for collision 



energies between 10 ^ and 10,000 cm ^ by applying both the close-coupling (CC) method 
and the coupled-states (CS) approximation. All th e CC and CS calculatio ns were performed 



using the nonreactive scattering code MOLSCAT (IHutson &: Green 



I994J ). In the quantum 



scattering calculations, the coupled-channel equations were integrated using the modified 
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log-derivative Airy propagator of I Alexander &: ManolopoulosI ( 119861 ) with a variable step 
size. The largest Legendre terms adopted in the potential expansion for H2 and CO were, 
respectively, 8 and 10. The numbers of Gauss integration points used in projecting angular 
components of the potential were 10, 11, and 12 for integration in 61, 62, and 0, respectively. 
The propagation was carried out to a maximum intermolecular separation of R = 100 A. At 
each energy, a sufficient number of total angular momentum partial waves was included to 
ensure convergence of the state-to-state cross sections to within 1%. The maximum value 
of the total angular momentum quantum number J employed in the calculations was 300. 
Since, the computation time for CC calculations scales as j^^^, the CS method is adopted 
for energies greater than 800 cm^^ which results in considerable time savings with a scaling 
proportional to Jmax; where jmax is the maximum si ze of the rotational basis. Further details 



about the scattering calculations can be found in 



Yang et al. 



feOOGaUbh. The adopted 



rotati onal constants for H2 and CO are 60.853 cm (iHuber fc Herzberg 



cm ^ (ILovas et al 



19791 ) and 1.9225 



for H2, and likewise j2 for CO. 
for the CO-H2 complex, V04, of 



19791), respectively. Hereafter, ji denotes the rotational quantum number 
n the present study, we adopt ed the four- dimensional PES 



Jankowski fc Szalewiczl (120051 ) 



Calculations of the collision energy dependence of state-to-state quenching cross 
sections for initial rotational states of CO, ^2=1; 2, ■ ■ -, 40, by collisions with both para- and 
ortho-H2 (i.e., ji = 0, 1, respectively) were performed. The cross sections were thermally 
averaged over the kinetic energy distribution to yield state-to-state rate coefficients from 
specific initial rotational states j2 as functions of the temperature T, 



k 



32^3: 



(8 \ 



(1) 



where crj^^f^ is the rotational transition cross section with j2 and J2 being respectively the 
initial and final rotational quantum number of CO, /i is the reduced mass of the CO-H2 
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complex, Ek the kinetic energy, and = {ksT) ^, where ks is the Boltzmann constant. 

4. RESULTS AND DISCUSSION 

All state-to-state cross sections and rate coefficients for quenching are available 
on the UGA Molecular Opacity Project website (www.physast.uga.edu/ugamop/). 



The rate coefficients are also available in t 
Molecular Database (LAMDA, ISchoier et al 



le format of the Leiden Atomic and 



20051 ) on our website and the LAMDA site 



(www.strw.leidenuniv.nl/~moldata/). Fig. [T] shows typical examples of state-to-state cross 
sections for j2 = 5 and 20 for collisions with para-H2. A number of resonances are evident 
between 10~^ and 10^ cm~^ which result in significant modulation in the rate coefficients. 
However, the magnitude of the re sonances ar e seen t o decrease with j2- Cross sections 



for some other j2 can be found in lYang et al. 



( l2006al Jbl). The agreement between the CC 



and CS calculations is shown to be excellent with a difference typically better than ~20%, 
justifying the adoption of the CS approximation at the higher energies. 

Quenching rate coefficients from initial rotational states, ^2=5, 10, 20, and 40, at 
temperatures ranging from 1 to 3000 K are shown in Figs. [2l-|9]for CO scattering with para- 
and ortho-H2. Unfortunately, we are unaware of any experimental rate coefficient data for 



these initial s tates. 
coefficients of Flower 



lerefo re, the current results are compared with the theoretical rate 



(l200l[ ). which were obtained over the limited te mperature range, 5 to 



Wernli et al. 



(12006r i 



400 K, with the Vgg PES. Similar results with V04 were obtained by 
for rotational levels of CO up to 5 and temperatures in the range 5—70 K. Wernli et al. 
also presented an analytic fit valid in the same temperature range. The rate coefficients 
calculated using their analytic relation are almost identical to the current results. 



For the quenching of j2=5 due to para-H2 collisions. Fig. [2] shows that for temperatures 
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between ~1 and 100 K, which is the van der Waals interaction-dominated regime, the rate 
coefficients exhibit an undulatory temperature dependence for Aj2 = —1, —2, and —3 
transitions due to the presence of resonances^ though the magnitude of the undulations 
decrease with increasing 72 • At temperatures above ~100 K, the rate coefficients generally 
increase with increasing temperature. Comparison with the rate coefficients of Flower 
shows that at temperatures higher than ~50 K, there is generally good agreement, except 
for J2 = and 2, where Flower's results get larger at temperatures above about 200 K. 
The results of Flower are smaller than the present rate coefficients for lower temperatures 
with the discrepancy increasing with decreasing temperature. Exceptions to this behavior 
are the transitions J2 = 5 — i- = 2 and 4, in which Flower's rate coefficients are generally 
larger than the present results. For quenching from initial levels j2 =10 and 20 given in 
Figs, m and SJ respectively, the state-to-state rate coefficients have similar structure to 
that shown for J2 = 5. The rate coefficients are dominated by small Aj2- For Aj2 = — 1 
transitions, a slight bump can be seen at temperatures between ~1 and 100 K for j2 = 10 
and 20 as shown in Figs. [3] and HI Figure [5] displays q uenching rate coefficients for j2 = 40. 



Flower! (120011 ) stopped at j2 = 29 for 



No previous results are available as the calculations of 
para-H2 collisions. 

For scattering by ortho-H2, the trends noted for para-H2 are also evident as displayed 

lose obtained for 



Flowei] fl200lr ) for 



in Figs. [6l|91 The quenching rate coefficients are of similar magnitude as t 
para-H2 and show very similar structure. The highest level considered by 
ortho-H2 collisions was j2 = 20. 

Comparing Figs. |2ll9l some general trends can be noted. With few exceptions, the 
state-to-state rate coefficients appear to attain constant values as they decrease to 1 K. 



^Additional typ ical examples of scattering resonances in cross section can be found in 



Yang et al 



fl2006a 



b|) along with rate coefficient comparisons for j2 < 4. 
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However, this is deceptive as an orbiting resonance is typically present in the cross 
sectio ns for collision ene rgies between ~0.01-1 cm~^ (see Fig. [T]) computed on the V04 



PES flYang et al 



2006al ). For temperatures less than shown, t he rate coefficients typically 



decrease (increase) prior to the Wigner regime (jWigneiill948l ) for intermediate (high) 



before eventually attaining finite limiting values, typically below ~1 mK. 

Disregarding the effect of the resonances, one would expect two possible trends in 
the relative ordering of the rate coefficient values a s a consequence of en ergy and angular 



momentum gaps. The exponential energy gap law fISteinfeld et al 



19911 ) predicts that the 



rate coefficients should decrease with increasing internal energy difference between the 
initial and final levels (i.e., with increasing |Aj2|)- However, CO is nearly homonuclear and 
rotational transitio ns in homonuclear molecules follow an even Aj2 propensity rule (e.g.. 



Shepler et al. 



20071 ). though the effect may diminish for large values of |Aj2|- Both effects 
appear to be important for CO-H2. For example, Aj2 = —2 is the dominant transition for 
quenching for temperatures greater than ~1000-3000 K. It is also the dominant transition 
for temperatures less than ~10 K, but only for j2 = 5 and 10. For other temperatures, 
Aj2 = — 1 transitions are predominant, while for some cases Aj2 = —3 is secondary, but 
with a rate coefficient larger than the Aj2 = —2 transitions. However, for Aj2 < —3, the 
exponential energy gap law appears to be fulfilled in nearly all considered cases. 



5. ASTROPHYSICAL IMPLICATIONS 

As discussed in the introduction, carbon monoxide plays an important role in a variety 
of astrophysical and atmospheric environments and has been detected in countless sources 
under a variety of excitation conditions. We discuss a few of these below focusing primarily 
on observed and predicted transitions from intermediate- and high-j2 in the FIR. 
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Rotational emission lines of CO have been observed in high-redshift [z) objects 
including the qu a sar S PSS J1148+5251 at z = 6.42 using the Very Large Array. 



Narayanan et al 



(I2OO8I ) modeled the emission from CO(f = 0, with j2 



10 with a 



non-LTE radiative transfer code and considered excitation due to H2 with rate coefficients 
taken from the LAMDA database. They found that the CO spectral energy distribution 
(SED) peaks at j2 = 8 at the beginning of the quasar phase, 2; ~ 7 — 8. Observations 
and modeling of the CO emission can serve to probe the hierarchical buildup of the host 
galaxy and the quasar phase and can be used to estimate the halo mass and quasar host 
morphology. 



Spaans fc MeijerinkI (120081 ) proposed that CO emission from j2 as high as 30 
could be observed from high- 2; massive accreting black holes using the Atacama Large 
Millimeter /submillimeter Array (ALMA). Such observations could provide a diagnostic of 
the radiation source: PDR or x-ray dominated region (XDR), for example. However, as the 
densities of such objects are typically 10^ cm~^, collisions with the dominant species, H2, 
will determine the CO rotational populations. In Fig. [101 the critical densities for CO due 
to para-H2 collisions are plotted (neglecting optical depth effects) where the critical density 



of rotational level j2 is defined as 

ncrij2) 



[e.g., 



Osterbrock fc Ferland 



(2) 



20061 ) and Aj^^ji are the spontaneous transition probabilities 



for dipole transitions. Rotational levels j2 ^ 10 are seen to be clearly out of thermal 
equilibrium for a gas density of 10^ cm"'^ requiring a non-LTE analysis . As the LAMDA 



database incorporates the CO-H2 rate coefficients of 



Wernli et al. 



(120061) and 



Flower 



(I2OO1I ) 



with extrapolations for para- and o rtho-H2 collisions above 72 



the predicted CO line intensities of 



Spaans fc MeijerinkI ( l2008l ) for emission from levels 



29 and 20, respectively. 



larger than J2 = 20 would likely be improved if the current rate coefficients are adopted. 
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Interestingly, their predicted CO SEDs peak at j2 = 10 and j2 = 25 for PDR and XDR 
environments, respectively. 



The ? o = 6 — »• 5 and j-? = 7 Q lines have been o 



NGC 253 ( 



Hai: 



Bradford et al. 



ev-Dunsheath et al. 



3 



2008 



Bradford et al. 



jsery ed in the starburst galaxy 



20031 ). Using a non-LTE analysis. 



20031 ) concluded that the cooling in the CO lines is so large that it must 



be balance d by a cosmic ray hea. ting rate ~800 times greater than that in the Milky Way. 



Previously, 



Krolik fc Lepp 



(Il989[ ) suggested that highly-excited j2 lines (e.g., j2 = 16 — 58) 
could be detectable from Seyfert galaxies giving diagnostics of the internal pressure of the 
dust-obscured torus. 

The Infrared Space Observatory studied a large variety of pre-main sequence objects 
(i.e., protostars) and detected numerous FIR CO rotational lines. For example, emission 



lines from j2 = 14 — 26 were observed in the Class object L1448-mm (INisini et al 



19991), 



j2 = 14-31 in the Class I object IC 139 6N, h = 14 



p = 14 — 20 in the Cla ss II object R CrA (jSaraceno et al. 



(ISpinoglio et al 



6 in the Class I object W28 A2, 



19991 ). and j2 = 14-25 in T Tauri 



19991 ). Non-LTE analyses of these pre-main sequence objects suggested 



multiple components with various excitation conditions. Future studies of these protostars, 
the infrared sources discussed above, and future FIR and submillimeter observations with 
Herschel and ALMA can benefit from the CO-II2 quenching rate coefficients presented here. 



6. CONCLUSIONS 



Rotational quenching of CO due to para- and ortho-H2 collisions has been studied 



using an explicit quantum-mechanical close-cou pling approach and the 



coup 



ed-states 



approximation on the potential surface, V04, of iJankowski fc Szalewiczl (120051 ). State-to- 
state quenching rate coefficients for initial rotational levels j2 =1? 2, ■ ■ -, 40 of CO are 
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obtained over a wide temperature range and available in tables formatted for astrophysical 
applications. Resonances result in undulations in the temperature dependence of the 
rate coefficients with the amplitudes of the undulations typically decreasing with j2- 
For temperatures less than ~50 K, the current state-to-state rotat ional quenchin g rate 



Flowei] (j200l|) obtained 



coefficients obtained with Vo4 are found to depart from the results of [ 
with the earlier Vgg PES. The discrepanices are likely related to the differences in the well 
depth and anisotropy of the two potentials. 
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Fig. 1. — Cross sections for the quenching of C0(j2) due to para-H2 coUisions as a func- 
tion of coUision energy obtained with the close-couphng method (hnes) and coupled-states 
approximation (symbols), (a) j2 = 5, (b) j2 = 20. 
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Fig. 2. — Rate coefficients for the quenching of CO(j2=5) by coUisions with para-H2 as 
functions of temper ature. L ines indicate present calculations on potential V04, symbols 
denote the results of 



Flowed (120011 ) on potential Vgg. Solid line and solid circles: ^2 = 5 — )■ 
j'2 = 0, dotted line and open circles: ^2 = 5 — )■ j'2 = 1, dashed line and solid squares: 
j2 = 5 — J- = 2, dot-dash line and solid squares: ^2 = 5 — )■ = 3, double-dot-dash line and 
solid triangles: j2 = 5 — )■ = 4. 
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Fig. 3. — Same as Fig. [21 except for j2 = 10 and as indicated. 
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Fig. 5. — Same as Fig. |2l except for j2 = 40 and as indicated. 
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Fig. 6.— 



Same as Fig. Ej except for ortho-H2. 



- 25 - 




Fig. 7.- 



Same as Fig. |3l except for ortho-H2. 
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Fig. 9.— 



Same as Fig. |5l except for ortho-H2. 
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Fig. 10. — Critical densities for CO{v = 0, j^) due to para-H2 collisions as a function of gas 
temperature T. 



